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摘 要 : 20 世纪 以 来 , 随 着 全 球 气候 变 暧 加 剧 , 冰 川 和 积 雪 普 讶 退缩 ,严重 影响 到 人 类 的 生存 和 社会 

经 济 的 可 持续 发 展 , 这 一 问题 在 我 国 西北 于 旱 区 的 博 格 达 峰 地 区 及 其 周边 地 区 尤为 突出 。 以 博 格 

达 峰 地 区 为 例 , 利 用 1990—2016 年 Landsat 5 5 Landsat 8 遥感 影像 ,对 比分 析 归 一 化 积 雪 指 数 
o (NDSI) 、 归 一 化 冰雪 指数 (NDS1) 、 归 一 化 主 成 分 雪 指 数 (NDPCS1) fa 2 tE He d 08 RE EC ( WET) Æ 
p 博 格 达 峰 地 区 监测 冰川 和 积 雪 的 能 力 , 同 时 结合 研究 区 周边 气温 、 降 水 数据 和 研究 区 地 形 数 据 , 探 
讨 博 格 达 峰 地 区 冰川 和 积 雪 面 积 交 化 与 区 域 地 形 、 气 候 间 的 响应 关系 。 结 果 表 明 :(1) eu] 
于 NDSII .NDSI $a NDPCSI 精度 值 更 高 ,可 以 替代 NDSI NDSII 监测 博 格 达 峰 地 区 冰川 和 积 雪 面积 。 
(2) 博 格 达 峰 地 区 冰川 和 积 雪 面积 呈 持 续 退 缩 的 趋势 。1990 一 2016 年 ， 052 22 aly: 
e 20.0796 , 且 年 退缩 率 不 断 增加 。(3) 高 程 .坡度 和 坡 向 对 冰川 和 积 雪 面积 变化 的 影响 较 显 著 , 山 地 
y 影 对 其 影响 较 弱 ,气温 的 升 高 是 冰雪 面积 减少 的 主要 因素 。 
m i $E ial: Landsat; 指数 ; 冰川 和 积 雪 变 化 ; 博 格 达 峰 
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CN 冰川 和 积 雪 是 反应 全 球 气 候 变 暖 的 指标 ,也 ， 研究 ,但 因 其 空间 分 辩 率 相对 较 低 (分 别 为 1.1 km 
— o0 是 干旱 地 区 人 类 赖 以 生存 和 发 展 的 重要 水 资源 ， 。 和 500 m) ,不 适用 于 区 域 研究 894) 。 近 年 来 ,Landsat 
对 区 域 生 态 环 境 、 人 民生 产生 活 和 社会 经 济 发 展 因 高 分 辨 率 (30 m) .全球 覆盖 和 丰富 的 数据 被 广泛 
起 到 至 关 重 要 的 作用 "i。 近 年 来 ,在 自然 和 社会 ” 用 于 区 域 冰川 和 积 雪 研究 '"。 在 冰川 和 积 雪 信 息 
— 经 济 因 素 双 重 驱 动 下 ,全 球 范围 内 冰川 和 积 雪 的 。” 的 提取 上 ,目前 主要 有 人 工 解 译 、 自 动 或 半自动 分 类 
G 数量 和 面积 也 呈现 出 不 同 程度 的 缩减 状态 = 。 ”以 及 基于 指数 的 提取 方法 5 。 人 工 解 译 和 自动 或 
长 期 的 冰川 和 积 雪 退 缩 可 能 会 使 区 域 水 资源 储量 。 半自动 分 类 虽然 精度 较 高 ,但 需要 丰富 的 判读 经 验 
降低 ， oe 这 必然 会 对 内 陆 ne 且 耗 时 耗 力 , 实 际 应 用 较为 困难 ”"” 。 基 
e de 9 影响。 因此 ， 指数 的 方法 通过 监测 阔 值 获取 冰川 和 积 雪 的 监测 
加 强 冰 川 和 积 雪 的 监测 E RE. 干 Rum ,由 于 其 快速 .简单 .准确 的 优点 ,在 研究 中 
旱 区 水 资源 变化 以 及 对 于 旱地 区 的 可 持续 发 展 至 。 得 到 了 广泛 的 应 用 。 如 BURNS 45", DURAN- 
关 重 要 。 ALARCON 等 使 用 归 一 化 积 雪 指 数 ( NDS1) 提取 
遥感 监测 是 研究 冰川 和 积 雪 的 普遍 方法 。NO- 秘鲁 科 迪 勒 拉 山脉 冰川 信息 ; SANKEY 等 ”、 
AA AVHRR(Advanced Very High Resolution Radiom- WALTERS 4&5! 使 用 NDSI Aan 言 息 IF 
eter) 和 MODIS( Moderate Resolution Imaging Spectro- 发 适合 的 积 雪 分 类 模型 和 降低 季 积 雪 融化 面积 
radiometer) 被 广泛 用 于 全 球 尺 度 上 冰川 和 积 雪 监 测 。 ”的 方法 ; 孙 永 猛 等 Pe de NDSI 提取 
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NDSI 的 准确 性 依赖 光谱 反射 率 ,受到 大 气 、 地 形 的 
影响 。KHOPKAR 45 R Jg IH — fe vk 5: 48 Zi 
CNDSIT) 提取 冰 川 信息 ,获取 冰川 时 空 分 布 变化 特 
征 及 原因 。SIBANDZE 57" 利用 主 成 分 分 析 
(PCA) 区 分 水 体 和 阴影 的 优越 性 开发 了 归 一 化 主 成 
分 雪 指 数 (NDPCSI) 。SATIR 将 织 帽 转换 湿度 指 
数 (WET) 用 于 冰雪 监测 ,发 现 比 NDSI 精度 更 高 ,但 
WET 难以 区 分 固态 和 液态 的 水 体 。 由 此 可 见 , 利 用 
遥感 览 测 冰川 和 积 雪 的 方法 有 很 多 ,但 不 同 研究 者 
基于 不 同 物理 指标 的 冰川 和 积 雪 监 测 方法 的 适用 性 
并 不 一 致 ,所 以 在 监测 区 域 ,对 冰川 和 积 雪 监测 方法 
进行 系统 验证 和 评估 ,确保 监测 结果 的 有 效 性 和 准 
确 性 是 必要 的 。 

博 格 达 峰 地 区 是 东 天 山 最 大 的 冰川 作用 区 ， 
受到 许多 学 者 的 关注 。 多 数 研究 者 提取 该 区 域 冰 
川 和 积 雪 信 息 时 主要 采用 比值 冰 值 法 和 NDSL。 如 
fa t^ 77 利用 比值 阔 值 法 提取 博 格 达 峰 冰 川 边 
界 ,获取 了 博 格 达 峰 冰川 持续 减少 的 信息 及 原因 。 
言 息 提 取 方 式 略 显 单薄 , 且 对 该 区 域 冰川 和 积 雪 
信息 提取 方法 准确 性 进行 评估 的 研究 不 足 。 分 析 
冰川 和 积 雪 覆盖 面积 变化 影响 因子 时 ,多 数 研究 
考虑 气候 因素 ,考虑 地 形 因素 影响 的 文章 较 少 。 
因此 ,本文 以 博 格 达 峰 地 区 为 研究 对 象 ,应 用 遥感 
和 地 理 信 息 系统 等 现代 手段 ,对 NDSI ,NDSII ND- 
PCSI WET 4 个 指数 的 冰川 和 积 雪 监测 能 力 进行 
了 评价 ,并 对 冰川 和 积 雪 面积 变化 及 其 原因 进行 
了 深入 研究 ,为 准确 提取 博 格 达 峰 地 区 冰川 和 积 
雪 信 息 以 及 周边 地 区 经 济 的 可 持续 发 展 提供 理论 
支撑 和 决策 参考 。 


1 研究 区 域 


博 格 达 峰 地 区 (88°03’ ~ 88°30’ E, 43°40! ~ 
43°55'N ,海拔 5 445 m) 位 于 新 疆 维吾尔 自治 区 息 康 
市 境内 ,是 天 山东 上段 博 格 达 山 的 主峰 (图 1)。 冬 季 
受 强 大 的 蒙古 冷 高 压 控制 ,寒冷 而 干燥 , 峰 顶 温度 常 
在 冰点 以 下 。 夏 季 高 空 盛行 西风 环流 ,是 典型 的 中 
纬 西 风 带 大 陆 型 气候 。 南 北 坡 年 平均 气温 0 7C 等温 
线 海拔 分 别 为 2 800 m 和 2 500 m^! 。 该 区 域 分 布 
有 113 条 现代 冰川 ,分 别 位 于 南 坡 的 古 班 博 格 达 河 
等 , 北 坡 的 三 工 河 等 河流 源头 。 
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1 博 格 达 峰 地 区 冰川 和 积 雪 地 理 分 布 
Fig 1 Geographic distribution diagram of Mount Bogda 
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glacier and snow 
2 数据 来 源 与 方法 


2.1 数据 来 源 

Landsat TM „Landsat ETM + 遥感 影像 数据 和 DEM 
数据 来 自 地 理 空 间 数 据 云 (http : // www. gscloud. 
cn/) 。 该 产品 经 过 了 大 气 校 正和 地 面 控制 点 几何 
校正 并 且 利 用 DEM 进行 地 形 校正 。 为 减少 不 同时 
期 冰雪 变化 所 造成 的 误差 ,选取 影像 时 间 间 隔 为 
10 a 左右 并 且 尽 可 能 选取 夏秋 季 获 取 的 少 云 或 者 无 
云 的 双 感 影像 ( 表 1)。 人 气象 数据 来 自 中 国 气象 数据 
网 (http: //data. ema. cn/site/index. html) ,选用 研究 
区 周边 的 吐鲁番 (站 点 编号 51573 ,89. 2°E ,42.9°N, 
海拔 345 m) 、 奇 台 ( 站 点 编号 51379, 89. 34^ E, 
44. O1°N ,海拔 793.5 m) 和 乌鲁木齐 (站 点 编 51463， 
87.37*E,43. 47? N ,海拔 917. 9 m)3 个 气象 站 点 近 
27 a 的 年 均 气 温和 年 均 降水 数据 。 
2.2 研究 方法 

本 研究 主要 分 3 个 步骤 进行 :(1) 使 用 大 气 校 
正和 几何 校正 后 的 2001 年 TM 遥感 影像 获取 NDS7、 

R1 影像 数据 信息 
Tab.1 Images used in this study 


PES ^ Ty — mm 0 PAE 

条 带 号 行 编号 ZE /m 
142 30 1990 - 10 -09 TM 0.44 30 
142 30 2000 -09 -02 TM 0.60 30 
142 30 2001 - 10 -23 TM 0.28 30 
142 30 2001 -09 -29 ETM + 0.06 30 
142 30 2010 -08 - 13 TM 0.03 30 
142 30 2016 -09 -30 ETM + 1.01 30 
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NDSII ,NDPCSI ET4 个 指数 ,并 利用 15 m 的 ETM + 
全 色 波 段 检测 精度 值 。(2) 使 用 精度 最 高 的 指数 监 
测 博 格 达 峰 地 区 冰川 和 积 雪 。(3) 分 析 冰 川 和 积 雪 
面积 变化 影响 因素 。 根 据 博 格 达 峰 地 区 的 实际 情况 
在 第 一 个 步骤 中 ,需要 使 用 高 分 辨 影像 检验 低 分 辨 
率 影 像 ,而 Landsat 8 与 Landsat 5 数据 在 2001 年 同 
时 存在 ,因此 本 文采 用 2001 年 数据 进行 精度 检验 
(图 2)。 

2.2.1 积 雪 和 冰川 监测 指数 ” 归 一 化 积 雪 指数 
(NDS7) ,通过 绿 波段 (0. 56 pum) 和 近 红 外 (1. 65 
um) 的 光谱 反射 率 差异 来 识别 冰川 和 积 雪 中 ; 归 一 
化 冰雪 指数 (NDS17) ,通过 红 波段 (0.63 ~0.69 jum) 
和 中 红外 波段 (1.55 ~1.75 pm) 光 谱 反 射 率 差异 来 
识别 冰川 和 积 雪 "i ; 归 一 化 主 成 分 雪 指 数 ( NDPC- 
SI) ,通过 对 landsat 波段 主 成 分 分 析 后 两 个 最 大 
PC 值 包 含 信息 的 差异 识别 积 雪 1 ; 继 帽 转换 湿度 
(WET) ,采用 BAIG 4&2?" ffj TM ETM + Bcd See $e 
换 方法 提取 WETness 数据 ,通过 结合 WETness 与 对 
应 波段 光谱 反射 率 识 别 冰 川 和 积 雪 ( 表 2) 。 

2.2.2 阅 值 与 精度 值 ” 通 过 精度 分 析 ,使 用 者 能 根 
据 分 类 结果 的 精度 ,正确 .有效 地 获取 分 类 结果 中 的 
信息 。 精 度 评价 主要 有 两 种 方案 :第 一 种 方案 是 
Landsat 产品 与 地 面 实 测 资料 进行 比较 ;第 二 种 方案 
是 与 更 高 分 辨 率 卫 星 影像 或 者 航 片 来 进行 比较 。 
在 缺少 站 点 的 山区 ,实测 数据 缺乏 且 成 本 过 高 , 且 获 
取 困 难 , 因 此 本 文采 用 第 二 种 方案 。 全 色 波 段 与 多 
光谱 波段 融合 后 空间 分 辩 率 为 15 m 的 ETM + 数据 
比分 辩 率 为 30 m 的 TM 数据 多 光谱 波段 优越 4 倍 ， 
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图 2 ”冰川 和 积 雪 检 监测 流程 
Fig2 Flow diagram of the methodology for glacier and 


snow monitoring 


R2 冰川 和 积 雪 监测 指数 


Tab.2 Glacier and snow monitoring indexes 


指数 计算 公式 
NDS1 9) (TM2 — TM5)/ ( TM2 + TM5) 
NDSII U9) ( TM3 - TM5 )/ ( TM3 + TM5) 
NDPCSIU9) — (1stPC -2ndPC)/(1stPC +2ndPC) 
Landsat 5 0. 150 9 x TM1 +0. 197 3 x TM2 +0. 329 9 x TM3 + 
WET?!) 0. 340 7 x TM4 -0. 711 2 x TM5 -0. 457 2 x TM7 
Landsat 8 0. 151 1 x OLD. +0. 197 3 x OLB +0. 328 2 x OLIA + 
WETO! 0. 340 7 x OLI50 -0.711 7 x OLI6 — 0. 455 9 x OLIT 


TE: TM 表示 Landsat 5 AIBOA; OLI 表示 Landsat 8 的 波段 
光谱 反射 率 ;PC 代表 主 成 分 值 


理论 上 可 以 用 来 评价 精度 。 为 了 方便 两 种 冰川 和 积 
雪 图 对 比 ,将 空间 分 辨 率 为 15 m 冰川 和 积 雪 图 重 采 
样 成 30 m 的 影像 。 利 用 混淆 矩阵 计算 总 精度 检验 
国 值 的 合理 性 。 


3 结果 与 分 析 


3.1 指数 阐 值 与 精度 值 确定 

通过 人 机 交互 式 技术 确定 4 个 指数 阔 值 的 取 值 
范围 ,为 减少 不 必要 的 计算 量 ,本 文 在 各 指数 阔 值 范 
围 内 ,以 0.05 为 间隔 ,获取 不 同 冰 值 的 冰川 和 积 雪 
图 ,统计 其 冰雪 和 象 元 个 数 ,分 别 与 Landsat ETM + vk 
川 和 积 雪 图 进行 对 比 确定 精度 最 高 的 闷 值 。 各 指数 
最 高 精度 的 混淆 矩阵 展示 在 表 3 中 。 

通过 总 精度 监测 , WET AE BE 0. 8 时 精度 最 高 
为 0. 94; NDPCSI X |j (EL 0. 55 时 精度 值 最 高 为 
0. 90; NDSII 1e [jfi 0. 45 时 精度 值 最 高 为 0. 92 ;ND- 
SI EBE 0. 4 时 精度 值 最 高 为 0. 91。WEkT 具有 最 
高 的 精度 值 (图 3)。 
3.2 博 格 达 峰 地 区 冰雪 面积 的 监测 

根据 在 博 格 达 峰 地 区 反复 验证 冰川 和 积 雪 覆 盖 
情况 ,得 到 精度 值 最 高 的 指数 WET, x WET 指数 应 
用 于 1990 年 .2000 年 .2010 年 .2016 年 的 遥感 数 
据 ,获取 冰川 和 积 雪 覆盖 图 , 计算 冰川 和 积 雪 面积 
( 表 4) 并 比较 面积 变化 情况 (图 4)。 

结果 表明 :1990 一 2016 年 间 冰 川 和 积 雪 面 积 整 
体 出 现 减少 趋势 , 且 速 度 明 显 加 快 。1990 年 9 月 ~ 
1999 年 9 月 退缩 率 约 7.32% ,2000 年 9 月 ~2009 
年 9 月 退缩 率 约 8.09% ,2010 4£9 H -2016 £9 H 
退缩 率 约 6. 17% ,1990 一 2016 年 冰川 和 积 雪 面积 退 
缩 率 约 20.07% 。 
3.3 影响 冰雪 面积 变化 的 因素 

根据 1990 一 2016 年 冰川 和 积 雪 变 化 情况 评估 
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Tabl.3 Confusion matrix between Landsat and Landsat-ETM + glacier and snow cover classifications 
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表 3 Landsat 与 Landsat-ETM + 冰川 和 积 雪 图 分 类 结果 混淆 矩阵 


指数 与 阔 值 混淆 和 矩阵 总 精度 Kappa 
Landsat ETM + 
NDSI =0.4 Landsat TM 冰川 和 积 雪 像 元 非 冰 川 和 积 雪 像 元 0.91 0. 64 
冰川 和 积 雪 像 元 182 309 132 474 
非 冰 川 和 积 雪 像 元 32 389 1 424 631 
NDSII =0.45 Landsat TM 冰川 和 积 雪 像 元 非 冰 川 和 积 雪 像 元 0.92 0. 68 
冰川 和 积 雪 像 元 188 588 119 377 
非 冰 川 和 积 雪 像 元 26 110 1 437 728 
NDPCSI =0. 55 Landsat TM 冰川 和 积 雪 像 元 非 冰 川 和 积 雪 像 元 0.90 0.56 
冰川 和 积 雪 像 元 140 248 101 654 
非 冰 川 和 积 雪 像 元 74 450 1 455 451 
WET =0.8 Landsat TM 冰川 和 积 雪 像 元 非 冰 川 和 积 雪 像 元 0.94 0. 69 
冰川 和 积 雪 像 元 131 209 54 329 
非 冰 川 和 积 雪 像 元 45 821 1 540 444 
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不 同 指数 的 冰川 和 积 雪 的 监测 结果 


Fig3 Snow and glacier detection results of different index 


影响 。 在 


研究 范围 选取 高 程 ,坡度 . 坡 向 .山地 阴影 数据 ,分 析 
RA 1990—2016 年 博 格 达 峰 地 区 冰川 和 积 雪 面积 变化 


Tab.4 Snow and glacier area variation of Mount Bogda 


from 1990—2016 


gg ee MED MPRA 
面积 / km E / km EAD 

1990 127.10 + 上 0.01 

1999 117.80+0.01 1990 一 1999 年 ”9.30 «0.02 7.32 

2009 108.27 x0.01 2000—2009 4F 9.53 +0.02 8.09 

2016 101.59 «0.003 2010—2016 ^F 6.68 +0.01 6.17 
1999—2016 4E 16.21 0.01 13.76 
1990—2016 年 25.51 x0.02 20.07 
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Fig 4 Continuous glacier and snow changes of Mount Bogda 
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十 格 达 峰 地 区 冰川 和 积 雪 随地 形 因素 的 变化 规律 。 
选取 研究 区 周边 3 个 气象 站 的 多 年 年 降水 和 年 均 气 
温 数据 ,分 析 博 格 达 峰 地 区 冰川 和 积 雪 随 气候 因素 
的 变化 规律 。 

3.3.1 高 程 的 影响 为 分 析 高 程 对 冰川 和 积 雪 覆 
盖 的 影响 ,本 文 以 25 m 为 间距 对 高 程 进行 了 分 类 ， 
共 202 类 。 通 过 对 研究 区 1990 年 9 H 2016 4E9 H 
的 冰川 和 积 雪 图 与 高 程 分 类 图 又 加 ,并 计算 各 高 程 
2016 年 相对 于 1990 年 冰川 和 积 雪 的 覆盖 率 变化 。 
冰川 和 积 雪 和 覆盖 率 定义 为 积 雪 和 冰川 面积 除 以 各 高 
程 研 究 区 面积 。 用 2016 年 减 去 1990 年 冰川 和 积 雪 
覆盖 率 代表 多 年 冰川 和 积 雪 变 化 率 , 变 化 率 大 于 0 
的 部 分 表示 冰川 和 积 雪 增加 ,反之 , 则 是 消融 ,计算 
结果 如 图 5 rn. 
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图 5 冰川 和 积 雪 和 覆盖 变化 率 与 高 程 的 关系 


Fig5 Variations of glacier and snow cover ratio with elevation 


图 5 显示 , 博 格 达 峰 地 区 冰川 和 积 雪 覆 盖 率 变 
化 随 海拔 总 体 呈 w 型 。2 500 m 之 前 和 4 500 mz 
后 趋 于 平缓 ,冰川 和 积 雪 和 覆盖 率 基 本 不 变 。2 500 ~ 
3 500 m 海拔 的 冰川 和 积 雪 开始 融化 ,消融 率 随 海 
拔 的 增 大 而 增 大 。3 500 ~3 800 m 海拔 ,冰川 和 积 
雪 消 融 变 慢 ,消融 率 随 海 拔 的 增 大 而 减少 。3 800 ~ 
4 000 m 海拔 冰川 和 积 雪 消融 速度 随 海拔 增加 而 加 
快 。4 000 ~4 500 m 冰川 和 积 雪 消 融 速 度 随 海拔 减 
Ro 
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满足 冰川 和 积 雪 的 温度 条 件 ,使 得 冰川 和 积 雪 覆盖 
的 决定 性 条 件 变 成 降水 ,而 新 疆 地 区 在 高 于 2 500 m 
海拔 的 区 域 降 水 量 随 高 程 增加 呈 近 似 二 次 曲线 类 
型 ,最 大 降水 高 度 带 在 山腰 (3 100 ~3 800 m) ,在 


5 500 m 左 右 出 现 第 二 极 大 值 带 一 -六 ,降水 抑制 了 
冰川 和 积 雪 的 消融 。 另 一 方面 , 随 着 海拔 高 度 的 增 
加 ,气温 下 降 , 水 汽 压 下 降 ,太阳 辐射 加 强 , 既 减少 了 
温度 引起 的 消融 又 增加 了 冰川 和 积 雪 的 升华 量 呈 ]。 
3.3.2 坡度 的 影响 ”为 了 分 析 坡 度 对 博 格 达 峰 地 
区 冰川 和 积 雪 冻 融 的 影响 ,本 文 以 1° 为 间距 对 坡度 
进行 了 分 类 , 共 90 类 。 与 前 文 高 程 对 冰川 和 积 雪 覆 
盖 的 影响 分 析 相 似 ,对 研究 区 1990 年 9 月 .2016 年 
9 月 的 冰川 和 积 雪 图 与 坡度 分 类 图 姜 加 得 到 冰川 和 


积 雪 覆盖 变化 率 ,结果 如 图 6 所 示 。 
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图 6 冰川 和 积 雪 覆盖 变化 率 与 坡度 的 关系 


Fig6 Variations of glacier and snow cover ratio with slope 


图 6 表明 , 博 格 达 峰 地 区 冰川 和 积 雪 变化 率 随 
WEE V 型 。 坡 度 低 于 40° 的 区 域 ,冰川 和 积 雪 和 覆 
盖 变 化 率 总 体 随 坡度 增 大 而 增 大 ,坡度 高 于 40° 的 
区 域 则 随 坡 度 增 大 而 减 小 ,最 后 趋 于 平稳 。 

冰川 和 积 雪 履 盖 率 随 坡度 呈现 非 线性 变化 ,在 
40° 前 变 率 增 大 ,40° 度 后 变 率 减 少 。 这 种 现象 一 方 
面 可 能 是 因为 不 同 坡度 条 件 下 因 正 午 太阳 高 度 导致 
太阳 辐射 和 气温 差异 ,引起 冰川 和 积 雪 升华 和 融化 
量 不 同 。 男 一 方面 ,在 坡度 低 于 40° 的 区 域 冰川 面 
积 随 坡度 增加 呈正 态 分 布 趋势 中 ,而 博 格 达 峰 地 区 
冰川 和 积 雪 主要 分 布 于 坡度 低 于 40° 区 域 。 

3.3.3 坡 向 的 影响 为 了 分 析 坡 向 对 博 格 达 峰 地 
区 冰川 和 积 雪 冻 融 的 影响 ,本 文 以 2° 为 间距 对 坡 癌 
进行 了 分 类 ,为 便于 统计 ,把 连续 的 坡 向 分 为 平地 
(一 1°)、 正 北 (337. 5° ~ 360°, 0° ~ 22. 5?) . 东北 
(22. 5° ~ 67. 5°)、 正 东 (67.5° ~ 112. 5°) ZR RH 
(112. 5° ~ 157. 5°) IE (157. 5° ~ 202. 5°) 西南 
(202. 5° ~ 247. 5?) | IE PĮ (247. 5° ~ 292. 5°) 西北 
(292.5° ~337.5°)9 类 。 对 1990 年 9 月 .2016 年 9 
月 的 冰川 和 积 雪 图 与 坡 向 分 类 图 琶 加 得 到 冰川 和 积 
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率 代表 多 年 冰川 和 积 雪 变化 ,结果 如 图 7 所 示 。 
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图 7 冰川 和 积 雪 覆盖 变化 率 与 坡 向 的 关系 


Fig7 Variations of glacier and snow cover ratio with aspect 


图 7 显示 , 博 格 达 峰 地 区 冰川 和 积 雪 和 履 盖 率 变 
化 总 体 呈 现 出 由 北 到 南 逐 渐 递 减 ,平地 区 域 的 变化 
率 较 低 。 从 北 坡 顺 时 针 至 西南 坡 , 变 化 幅度 不 断 减 
小 ,由 西南 坡 至 北 坡 逐 渐 增 大 , 正 北 方向 变化 率 最 
高 ,西南 方向 的 变化 率 最 低 。 
韦 格 达 峰 地 区 冰川 和 积 雪 履 盖 率 变化 随 坡 向 
(由 西北 向 西南 ) 呈现 逐渐 变 小 的 趋势 。 这 种 现象 
可 能 是 坡 向 影响 水 汽 运动 方向 和 太阳 光照 引起 温度 
差异 造成 ”] 。 博 格 达 峰 地 区 位 于 盛行 西风 带 , 水 汽 
主要 来 源 于 偏 西方 向 ,西南 光照 时 间 和 太阳 高 度 均 
大 于 东北 向 。 水 汽 和 太阳 光照 引起 的 冰川 和 积 雪 分 
布 不 均 。 
3.3.4 山地 阴影 的 影响 “为 了 分 析 山 地 阴影 对 博 
格 达 峰 地 区 冰川 和 积 雪 冻 融 的 影响 ,本文 以 2 为 间 
距 对 山地 阴影 灰 度 值 进行 了 分 类 , 共 127 类 。 对 研 
究 区 1990 年 9 月 .2016 年 9 月 的 冰川 和 积 雪 图 与 山 
地 阴影 分 类 图 释 加 得 到 冰川 和 积 雪 和 覆盖 率 ,用 2016 
年 减 去 1990 年 冰川 和 积 雪 覆 盖 率 代表 多 年 冰川 和 
积 雪 变 化 ,结果 如 图 8 所 示 。 

图 8 表明 , 博 格 达 峰 地 区 冰川 和 积 雪 禾 盖 变化 
率 随 坡度 上 下 波动 ,整体 上 随 山 地 阴影 灰 度 值 增加 
而 消融 速度 变 绥 。 但 山地 阴影 灰 度 值 100 左右 冰川 
和 积 雪 消融 最 快 ,山地 阴影 灰 度 值 230 左右 冰川 和 
积 雪 消 融 缓慢 ,规律 性 不 明显 。 因 此 ,我 们 可 以 推测 
山地 阴影 与 冰川 和 积 雪 覆盖 率 变 化 关系 不 显著 ,对 
冰川 和 积 雪 冻 融 的 影响 较 小 。 
3.3.5 气候 的 影响 为 了 分 析 气 候 对 博 格 达 峰 地 
区 冰川 和 积 雪 的 影响 ,本 文 结合 博 格 达 峰 地 区 边缘 
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图 8 冰川 和 积 雪 覆盖 变化 率 与 山地 阴影 的 关系 


Fig8 Variations of glacier and snow cover ratio with hillshade 


的 奇 台 ,吐鲁番 、 乌 鲁 木 齐 3 个 气象 站 的 降水 年 距 平 
和 年 均 温 变 化 数据 (http: // data. cma. en/) ,推测 出 
村 格 达 峰 地 区 降水 年 距 平 和 年 均 气 温 变化 量 ( 图 
9)。 张 明 军 等 ”认为 气温 每 升 高 1 必需 要 降水 量 
增加 47% 或 56% 才 能 弥补 由 升温 引起 的 冰川 消融 。 
本 文选 取 3 个 站 点 的 平均 升温 率 大 约 在 0.30 "C. - 
(10 a) 左右 ， 而 降水 增加 幅度 只 有 13 mm - 
(10 a) 左右 ,降水 量 增加 远 不 足以 弥补 升温 引起 
的 消融 。 何 毅 等 5 利用 Morlet 小 波 分 析 研 究 表明 
东 天 山 博 格 达 峰 地 区 在 1972—2013 年 处 于 气温 偏 
高 降水 偏 少 期 。 因 此 ,可 以 推断 气温 对 博 格 达 峰 地 
区 冰川 和 积 雪 面 积 变 化 影响 较为 显著 。 


4 讨论 


本 研究 提出 基于 指数 的 冰雪 监测 方法 ,通过 对 
志 格 达 峰 地 区 冰川 和 积 雪 情况 监测 反复 验证 NST , 
NDSII NDPCSI WET 的 适用 性 。 本 研究 中 ,WET 指 
数 被 证 明 是 最 准确 的 指标 ,但 WET 和 NDSI 的 准确 
性 差异 极 小 ,很 容易 受 地 物 分 布 状况 、 云 量 等 因素 的 
影响 ,导致 准确 性 发 生变 化 。SIBANDZE 等 "在 研 
究 祁连山 的 积 雪 时 发 现 VDPCS1 监测 积 雪 比 NDSI 
有 更 好 的 精度 。 但 本 文 NDPCSI 监测 最 高 精度 是 
0.9, 低 于 NDSI, NDSII WET 的 最 高 精度 。 在 SAT- 
IR 的 研究 中 也 出 现 了 类 似 的 情况 ,这 可 能 是 受 数 
据 类 型 或 冰雪 的 差异 影响 。 

本 文 使 用 WET 指数 监测 博 格 达 峰 地 区 冰川 和 
积 雪 变 化 时 发 现 1990—1999 年 .2000 一 2016 年 的 
退缩 率 与 何 毅 等 ”5 的 研究 结果 基本 一 致 ,表明 使 
用 WET 指数 监测 博 格 达 峰 地 区 冰川 和 积 雪 是 可 行 
的 。 而 差异 的 出 现 可 能 是 受 研究 区 域 范 围 、 冰川 数 


周 远 刚 等 : 博 格 达 峰 地 区 冰川 和 积 雪 变化 遥感 监测 及 影响 因素 分 析 


—o— FHC 


一 人 -一 降水 年 距 平 /mm 


要 包括 地 形 和 气候 两 个 方面 。 张 明 军 等 ” 认为 气 
候 是 影响 博 格 达 峰 地 区 冰川 和 积 雪 面积 变化 的 主要 
因素 ,但 本 文 在 分 析 博 格 达 峰 地 区 附近 的 奇 台 气象 
站 的 气温 .降水 变化 时 发 现 1990 一 1999 年 的 平均 年 
均 温 为 5.63 C ,平均 年 降水 量 为 221. 1 mm,1999— 
2013 年 的 平均 年 均 温 为 5.57 7C ,平均 年 降水 量 为 
196.85 mm。 冰 川 和 积 雪 面积 退缩 与 温度 和 降水 变 
化 线性 关系 不 明显 。 因 此 地 形 因素 在 博 格 达 峰 地 区 
冰川 和 积 雪 面 积 变化 中 有 着 重要 作用 。 尽 管 气温 是 
冰川 和 积 雪 冻 融 的 关键 因素 ,但 坡度 . 坡 向 等 也 有 重 
要 的 影响 |。 气 候 因素 与 地 形 因 素 结合 ,更 多 的 大 
阳 辐 射 和 更 少 的 降水 量 能 引起 更 多 的 冰川 和 积 雪 融 
化 和 升华 。 而 山地 阴影 并 没有 像 预 期 的 那样 影响 冰 
川 和 积 雪 冻 融 ,这 意味 着 太阳 辐射 并 不 是 影响 博 格 
达 峰 地 区 冰川 和 积 雪 的 重要 因素 。 


5 结论 


本 文 利 用 博 格 达 峰 地 区 1990 一 2016 年 的 Land- 
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图 9 年 均 气 温 与 年 降水 中 平 变化 
Fig9 Variation of average annual temperature and 


precipitation anomaly 


目 和 分 布 区域 的 差异 等 的 影响 。 对 比 近 几 年 我 国 西 
部 冰川 变化 研究 成 果 : 祁 连 山 中 段 冰川 年 退缩 率 为 
0.46% :al; 阿 尔 金 山地 区 冰川 年 退缩 率 为 
0.42% -a- UP ;喜马拉雅 山东 段 治 扎 冰 川 年 退缩 
率 为 0.60% - a !U* ;克拉 昆仑 山 克 勤 青 河流 域 冰 
川 年 退缩 率 为 0.22% . a CUE, GER, TES 
球 变 暖 的 大 背景 下 , 博 格 达 峰 地 区 冰川 和 积 雪 变化 
趋势 与 中 国 西 部 其 他 冰川 所 反映 出 来 的 大 趋势 相 吻 
合 ,但 退缩 速率 上 存在 差异 ,这 可 能 和 冰川 面积 、 盛 
行 风 影响 等 有 关 。 

韦 格 达 峰 地 区 冰川 和 积 雪 面 积 变化 影响 因素 主 


sat TM ETM + 


数据 分 析 了 NDSI、NDSIT、NDPCST、 


WET 4 个 指数 的 适用 性 ,以 及 地 形 和 气候 因素 对 冰 
川 和 积 雪 冻 融 的 影响 ,得 到 了 以 下 主要 结论 : 

(1) Wk7T 能 很 好 的 蔡 代 NDSI 和 NDSI 精确 监 
测 博 格 达 峰 地 区 冰川 和 积 雪 。 在 博 格 达 峰 地 区 反复 
验证 NDSI, NDSII , NDPCSI WET 的 适用 性 ,4 个 指 
数 的 精度 表现 为 :WET > NDSII > NDSI > NDPCSI , 

(2) 1990—2016 年 , 博 格 达 峰 地 区 冰川 和 积 雪 
面积 持续 减少 , 且 速 度 呈 增加 趋势 。1990 一 1999 年 
冰川 和 积 雪 面积 退缩 率 约 7.32% ;2000 一 2009 年 冰 
川 和 积 雪 面积 退缩 率 约 8. 09% ;2010 一 2016 年 冰川 
和 积 雪 面积 退缩 率 约 6.17% 。 

(3) 冰 川 和 积 雪 气 候 和 地 形 因素 呈现 出 不 同 的 
特点 。 地 形 要 素 中 ,冰川 和 积 雪 覆 盖 变 化 率 总 体 随 
海拔 增 大 而 增 大 , 随 坡度 增 大 呈现 V 型 分 布 , 随 坡 
向 变化 由 西南 向 东北 方向 减 小 。 气 候 要 素 中 ,气温 
是 冰川 和 积 雪 履 盖 变 化 的 主要 因素 ,降水 影响 较 弱 。 
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Remote sensing monitoring of the change of glacier and snow cover 
and its influencing factors in Mount Bogda 


ZHOU Yuan-gang ^, ZHAO Rui-feng ^, Zhang Li-hua'^, ZHAO Min'? 
(1 College of Geography and Environment Science , Northwest Normal University , Lanzhou 730070 , Gansu , China ; 
2 Gansu Engineering Research Center of Land Utilization and Comprehension Consolidation , Lanzhou 730070 , Gansu , China ) 


Abstract: Glaciers and snow cover are indicators of global warming and important water resources for human sur- 
vival and development in arid regions. Since 1900s ,the glacier and snow area has generally been shrinking as global 
warming intensifies , which seriously affects people 's life , social and economic activities in the related areas. Particu- 
larly ,in the surrounding areas of Mount Bogda which located in the arid area of northwest China, this kind of effects 
is much more significant. By strengthening the monitoring the change of glacier and snow cover, we can not only un- 
derstand the characteristics and trends of glacier and snow change , but also reveal the response of glacier and snow 
to climate change , and thus to possibly predict future changes of glacier and snow. We chose Mount Bogda at Fukang 
County , Xinjiang , China as the research area and obtained remote sensing images of Landsat 5 and Landsat 8, DEM 
data and surrounding temperature and precipitation data from 1990 to 2016. Firstly , we employed the satellite images 
of Landsat 5 and Landsat 8 to calculate value of Normalized Difference Snow Index ( NDST) , Normalized Difference 
Snow-Glacier Index ( NDSII) , Normalized Difference Principle Component Snow Index ( NDPCSI) and Tasseled 
Cap Wetness transformation ( WET) and to analyze the respective ability to detect glacier and snow. Secondly , we 
used WET with high accuracy to obtain glacier and snow cover data in 1990 ,2000 ,2010 and 2016 respectively and 
analyzed the changes of glacier and snow cover. Finally , combining with the temperature data and the terrain data of 
the study area, we explored the relationship between the change of snow and glacier area of Mount Bogda and the 
terrain and climate in this area. The results show as follows; (1) The WET index has a higher overall accuracy than 
the NDSI index , which can replace NDSI and NDSII to detect and monitor the glacier and snow area changes in 
Mount Bogda. (2) The area of snow and glacier continued to shrink in Mount Bogda. From 1990 to 2016, the snow 
and glacier area was reduced by 20.07% ,and the annual reduction rate continues to increase. (3) The influence 
of temperature and slope on the change of snow and glacier area is strong, while the influence of mountain shadow 
and altitude is weak. And the rising temperatures is the main factor of the decreased area of snow and glacier. The 
research results could provide theoretical support and decision reference for accurate extraction of glacier and snow 
cover information in Mount Bogda area and for the sustainable economy development in the surrounding areas. 


Key words: Landsat; index; variations of glacier and snow; Mount Bogda 


